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We explore the photon statistics of light emitted from a system comprising a single four-level atom 
strongly coupled to a high-finesse optical cavity mode which is driven by a coherent laser field. In 
the weak driving regime this system is found to exhibit a photon blockade effect. For intermediate 
driving strengths we find a sudden change in the photon statistics of the light emitted from the 
cavity. Photon antibunching switches to photon bunching over a very narrow range of intracavity 
photon number. It is proven that this sudden change in photon statistics occurs due to the existence 
of robust quantum interference of transitions between the dressed states of the atom-cavity system. 
Furthermore, it is shown that the strong photon bunching is a nonclassical effect for certain values 
of driving field strength, violating classical inequalities for field correlations. 

PACS numbers: 42.50.-p, 32.80.-t, 42.65.-k 



I. INTRODUCTION 



A promising avenue of research in the attempt to build 
quantum computers uses large optical nonlinearities to 
create essential elements for quantum computation, such 
as quantum gates pL Obtaining large, noiseless non- 
linearities is of course a fundamental goal in the field 
of nonlinear optics, and technological advances in recent 
years, for example in cavity quantum electrodynamics 
(CQED) [|| and the ability to access and study strongly 
coupled quantum systems, offer an exciting new era in 
this field. 

Attempts to achieve large optical nonlinearities are 
generally plagued by a seemingly necessary trade-off be- 
tween large absorption and the strength of nonlincar- 
ity: to obtain a strong nonlinearity one has to drive the 
atomic system close to resonance, triggering large absorp- 
tion and therefore a strong source of noise. However, it 
is possible to reduce atomic absorption (and hence spon- 
taneous emission) by utilizing quantum coherence effects 
in multilevel atoms. In particular, electromagnetically in- 
duced transparency (EIT) |3) was employed by Schmidt 
and Imamoglu to devise a scheme involving four-level 
atoms which produces a large Kerr nonlinearity with vir- 
tually no noise. It was then shown by Imamoglu et al. ^ 
that if such a strong optical Kerr nonlinearity is imple- 
mented in a CQED setting, then it is possible to realize 
the effect of photon blockade, in which the atom-cavity 
system effectively acts as a turnstile device for single pho- 
tons. The physical explanation of this effect is simple. 
Only the transition between the ground and first excited 
state of the nonlinear atom-cavity system is resonant 
with the driving field. Higher states are detuned from 
resonance by an amount proportional to the strength of 
nonlinearity. Tian and Carmichael |6j have also predicted 
such an effect with a single two-level atom strongly cou- 
pled to the cavity mode. 
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The proposal of Schmidt and Imamoglu [|f|, although 
very appealing in its use of EIT to reduce decoherence, 
meets obstacles in the attempt for implementation in the 
many-atom CQED environment. It was shown that this 
particular realization of the photon blockade system is 
not ideal for demonstrating photon blockade, because of 
the strong linear dispersion of the medium fil g[. In an 
attempt to avoid this difficulty, Rebic et al. B proposed 
a model in which a single four level atom is trapped in a 
high-finesse microcavity. It was shown that this system 
(called the single-atom EIT-Kerr system) can act as a 
near-ideal Kerr optical nonlinearity. In such a strongly 
coupled system, the excitations can be labeled as polari- 
tons, which are defined as mixtures of atom/cavity mode 
excitations. For weak to moderate driving, the EIT-Kerr 
system can be approximated by a two-state system, cor- 
responding to the two lowest lying polariton eigenstates. 
The transitions to the next set of excited states (the sec- 
ond manifold) , obtained by introducing a second photon 
from the driving field into the cavity, are highly detuned 
from the bare-cavity resonance, and therefore cannot be 
excited. Hence, in effect, the weakly driven single-atom 
EIT-Kerr system acts as an ideal photon blockade device. 

In a further work Gq], a Hamiltonian for the effective 
two-level system was derived using a polariton basis, and 
it was shown that such a Hamiltonian can describe the 
spectral properties and energy splittings (dynamic Stark 
effect) seen in the full model. Furthermore, to explain 
the properties of the second order correlation function, it 
was shown that the effective two-level system has to be 
generalized to include two extra states in the first excita- 
tion manifold [ p"l| ] . If more than one atom is introduced 
the situation becomes more complicated, since each addi- 
tional atom introduces extra energy levels into the second 
manifold, some of which are connected to the lower levels 
via transitions which are resonant with the cavity mode. 
However, it was shown by Werner and Imamoglu |T^ ] 
that the introduction of an additional atomic detuning 
can be used to solve this problem (see also the work of 
Greentree et al. JT^j). 

In |llj , a comparison of the EIT-Kerr scheme and the 
extended Jaynes-Cummings scheme (i.e. the standard 
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Jaynes-Cummings model |1J] with the effects of losses 
included) was made in terms of their effectiveness in pro- 
ducing photon blockade. The EIT-Kerr scheme exhibited 
a value of </ 2 ) (0) (the second order correlation function at 
zero time delay) of that of the corresponding Jaynes- 
Cummings scheme. This was attributed to destructive 
quantum interference between certain transition ampli- 
tudes. In particular, with a suitable choice of parame- 
ters, probability amplitudes for transitions from the first 
excited state (in the first manifold) to the two second 
manifold eigenstates closest to the cavity resonance can- 
cel each other. This leads to enhanced antibunching of 
the photons leaving the cavity so that g^(0) ~ 0. 

We have remarked that the validity of these results de- 
pends on the assumption of weak driving. It is therefore 
of interest to explore how robust the photon blockade 
is as the driving field strength is increased. We have 
performed an analysis of the driving field dependence in 
the EIT-Kerr system, and briefly compared the results 
with those obtained for the extended Jaynes-Cummings 
model. A surprising feature is found in the EIT-Kerr 
system, where a sudden change in photon statistics (i.e. 
of <7^(0)) of light emerging from the cavity occurs at 
intermediate driving strengths. If plotted against the 
mean intracavity photon number, g^(0) is seen to un- 
dergo a sudden transition, switching rapidly between an- 
tibunching and strong bunching. The exact position of 
the threshold depends on the characteristic system pa- 
rameters, namely the atom-field interaction strength and 
the (EIT) coupling field Rabi frequency. 

In Section II we present our model and the methods 
of solution employed. Section III contains the results of 
numerical simulations for the photon statistics, showing 
the surprising behaviour of the second order correlations 
in the regime of intermediate driving. Section IV con- 
tains a detailed explanation for the results presented in 
the preceding Section, obtained by analyzing the density 
matrix expressed in an appropriate basis. In Section V 
we present analytical calculations based on insights ob- 
tained in Section IV. 



II. MODEL 

A. Hamiltonian 

The driven atom-cavity configuration and the atomic 
energy level scheme are shown schematically in Fig. [|. 
The Hamiltonian describing the system in the rotating 
wave and electric dipole approximations is H = Ho + Hd, 
where, in the interaction picture, 

Ho = h5 ct 2 2 + fi-A (T44 + ihgi (aVi 2 — 021a) 

+iMl c (023 - 0-32) + ihg 2 (a t fr 34 - (743a), (la) 
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In these equations, OYj are atomic raising and lowering 
operators (for i ^ j), and atomic energy level population 



Driving £ p 



-H i-'i I 



9 (2) (r) 



('<) 




FIG. 1: (a) Envisaged experimental setup. A single four-level 
atom is trapped in an optical cavity with the decay constant 
k. The cavity is driven by a coherent laser field which couples 
to the cavity mode with strength £ p . An additional laser with 
Rabi frequency £l c directly couples to the atomic transition. 
(6) Atomic energy level scheme. The cavity mode drives tran- 
sitions |1) — > 1 2) and |3) — > j4), with the respective coupling 
strengths g\ and <?2- The transition j2) <-> |3) is coupled by a 
classical field of frequency uj c and Rabi frequency Q c . Spon- 
taneous emission rates are denoted by 7,- . Detunings 5 and A 
are positive in the shown configuration. 



operators (for i = j); a and are the cavity field annihi- 
lation and creation operators, respectively. Detunings 5 
and A are defined from the relevant atomic energy levels; 
gi.2 are atom-field coupling constants for the transitions 
|1) — ► 1 2) and |3) — > |4), and il c is the Rabi frequency of 
a coherent field coupling the transition |2) — > |3). The 
cavity driving field is characterized through the param- 
eter S p , related to the power output of the driving laser 



V via 



VnT 2 



(2) 



where T is the cavity mirror transmission coefficient, k 
is the cavity decay rate, and ui cav is the cavity mode fre- 
quency. Damping due to cavity decay and spontaneous 
emission is introduced below in the context of the quan- 
tum trajectory approach fl5f| . 

Four atomic levels plus the cavity mode span a Hilbert 
space of dimension 4 x N, where N denotes the trunca- 
tion order in a Fock state expansion of the cavity field 
subspace. In the absence of driving, Hamiltonian 7io, 
given by Eq. (|Ta|) , takes a block-diagonal form, with TV 
blocks on the main diagonal. Each block represents a 
manifold of eigenstates associated with the appropriate 
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term in the Fock expansion. The ground, first and sec- 
ond manifolds have been analyzed from the viewpoint of 
photon blockade in Refs ||, [l^, The general struc- 
ture of the dressed states in an arbitrary n th manifold 
has been discussed in Ref. fLOf . 

Addition of the driving term (^) complicates the anal- 
ysis, since the Hamiltonian matrix does not retain the 
block-diagonal form. It is possible, however, to re- 
express the driving Hamiltonian in terms of polariton 
operators and effective Rabi frequencies for transitions 
between dressed states. This has been done in pX| , and 
in the remainder of this article we will draw on these 
results. 

Our analysis of the driven 'atom-cavity molecule' will 
proceed in two complementary directions. First, using 
the method of quantum trajectories fi"5|| , we obtain nu- 
merical data. Then, using the formalism of Ref. , we 
construct an effective Hamiltonian in the polariton basis, 
which encapsulates the physics sufficiently to explain the 
numerical results. 



B. Damping: Quantum Trajectories 

In the quantum trajectories approach, damping enters 
the model through collapse operators, with each of these 
corresponding to one decay channel. In the case under 
consideration we require the following four collapse op- 
erators, 



C3 = \/73 cr 34 , Ci^\fna, 



(3) 



where 7^ denote spontaneous emission rates into each 
of the decay channels, and k is the cavity field decay 
rate. The effective Hamiltonian used in the trajectories 
approach is non-Hermitian and takes the form 



4 
k=l 



(4) 



with H given by ([l]). 

In deciding on the truncation for the cavity mode 
Hilbert space, special care must be taken to include a 
sufficient number of states to capture all the relevant dy- 
namics. If an empty cavity would be driven by an exter- 
nal coherent field corresponding to the parameter £ p , the 
amplitude of the intracavity coherent field would be 



iSp/ k , 



(5) 



and the expected mean intracavity photon number (n) = 
\a\ 2 . Our simulations include driving strengths of up to 
£ p = 3k, so the truncation is set at N = 40. The inset 
in Fig. g shows that the actual mean intracavity photon 
number stays well below its empty cavity counterpart 
for the given range of driving, thus justifying the chosen 
truncation. 



III. NUMERICAL SIMULATIONS OF PHOTON 
STATISTICS 

In this section we present the results of numerical sim- 
ulations using the method of quantum trajectories |fl5| , 
with results averaged over 10000 trajectories. Values of 
parameters chosen for the simulations are consistent with 
recent experiments [^6|, so the scheme presented in this 
paper should be experimentally viable. 

We evaluate the second-order correlation function 
<T 2 )(0) as a function of driving strength £ p . It was estab- 
lished earlier ^) that this function is a good measure 
of photon blockade; perfect photon blockade yields per- 
fectly antibunched photons. The steady-state second- 
order correlation function is given by 



(2) = ( t(t) t(t + T )a(t + T)o(t)) 

9 [> t-.oo(ot(i)a(t))(at(i + T)o(t + r)) 



(6) 



Perfect antibunching/photon blockade corresponds to 
g (2) {T = 0) = 0; for a coherent field gW(r = 0) = 1, 
and for a bunched field g (2 >{T = 0) > 1 @. 

Fig. I shows (t = 0) as a function of the (scaled) 
driving parameter for a single-atom EIT-Kerr system, 
compared with the same quantity calculated for the 
extended Jaynes-Cummings model. The extended JC 
model is driven on the lower Rabi resonance, as envis- 
aged by Tian and Carmichael ||. In the weak driving 
regime, simulation confirms earlier results - stronger pho- 
ton blockade in the EIT-Kerr system. As the driving in- 
creases, the extended Jaynes-Cummings model gradually 
saturates, and the field correlation asymptotically tends 
to one. Naively, one would expect qualitatively simi- 
lar behavior in the EIT-Kerr model. Our simulation, 
however, shows a vastly different result. Over a narrow 
range of driving, the statistics of the field changes from 
strongly antibunched to strongly bunched, and the coher- 
ent state value <?( 2 ' (0) = 1 is approached asymptotically 
from above as £ p is increased further. 

Note that the quantity being increased here is the 
scaled driving parameter, so one may argue that it is 
the ratio that determines this behavior, i.e. we may ei- 
ther increase the driving strength or decrease the cavity 
decay rate to obtain the same result. However, to clarify 
this issue as related to the photon statistics, in Fig. [| we 
show a parametric plot of second-order correlation func- 
tion against the expectation value of intracavity photon 
number, (a^a) ss . Both quantities are now unsealed by 
any external parameter, and the sudden nature of the 
switch becomes even more obvious, showing a phase- 
transition-like behaviour. This can be related back to 
Fig. |^. In particular, the transition from antibunching 
to bunching happens over a range of driving for which 
the intracavity photon number stays practically constant 
(see inset of Fig. ||), leading to the suddenness of the 
transition seen in Fig. |^, and the concentration of points 
around the threshold region. To further emphasize this, 
111 Fig. I (b) we plot the ' pure' second-order correlation 
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FIG. 2: Second order correlations at zero time-delay against 
the (scaled) driving parameter. The inset shows mean intra- 
cavity photon number in the steady state. The dotted line 
denotes the extended Jaynes-Cummings scheme for which 
7 = 0.1k, g — 6k. Other curves show the single-atom EIT- 
Kerr model with — 0.1k, gi = 6k, 5 = —0.2k, A = 0.1k. 
The solid line represents the case Q c — 2k; dot-dashed line 
Q. c = 6k; dashed line fi c = 18k. All curves are obtained by 
averaging over 10 4 trajectories. 

function (aWaa) ss which also exhibits a threshold- like 
behaviour. 

One other feature of the numerical results is noted. 
The bunching-antibunching transition is sharper and the 
subsequent bunching stronger for smaller Q c . In fact, it 
is the increasing ratio gi/ft c that really matters. The 
sharper transitions also occur at smaller values of (a'a) ss . 
For a decreasing ratio gi/Q c , </ 2 )(0) approaches the be- 
haviour of the extended JC model. 



IV. DENSITY MATRIX TREATMENT 

We proceed to determine which eigenstates of the 
strongly coupled quantum system contribute significantly 
to its dynamics. The total density matrix of the system 
in the steady state can be written in most general terms 
as 

p = J2cab\a){b\, (7) 

a, 6 

where a, b belong to the set of all possible states of the 
system in an arbitrary basis. The natural basis for the 
simulation is the one of bare states. The cavity mode 
subspace is truncated at 40, making the size of p equal 
to 160 x 160. Let T be the transformation that diago- 
nalizes Hamiltonian (|l]), i.e. the full Hamiltonian, with 
driving included. The density matrix can be transformed 
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FIG. 3: Second order correlations at zero time-delay against 
the mean steady-state intracavity photon number for a 
single-atom EIT-Kerr system. Parameters are the same as 
in Fig. ^| Figure (a) shows normalized correlation function 
with main figure showing the range where the sudden change 
can be seen in detail, while the inset shows the whole range. 
Figure (6) shows the unnormalized second order correlations. 



into a new basis as pdiag = TpT -1 , and we can look for 
the nonzero elements of this matrix. Diagonal elements 
of the matrix correspond to populations of the dressed 
states, while off-diagonal elements correspond to coher- 
ences between the dressed states. The results depend on 
the size of £ p , i.e. we expect the number of non-zero el- 
ements to increase as £ p is increased. For a large part of 
the £ p range considered, however, the density matrix is 
dominated by elements corresponding to two particular 
states. 

Figure |] shows the nonzero matrix elements of den- 
sity matrix pdiag- Diagonal states are sorted in increas- 
ing order; the state with largest negative detuning is at 
(h j) = (1) 1)) while the state with largest positive de- 
tuning is at (i, j) = (160, 160). States with the smallest 
detuning (i.e. closest to the cavity resonance) are at the 
centre of matrix, at positions 80 and 81 along the main di- 
agonal. Two diagonal elements dominate the matrix, and 
we identify these as being the populations of the Stark- 
split states \ip±) (see Appendix |A| and [[l0|). Stark-split 
states are therefore found at the positions 80 and 81 on 
the main diagonal of the density matrix. At positions 





FIG. 4: Relevant density matrix elements at the different values of driving. Elements of the central submatrix (dimension 
21 x 21) of the total density matrix (160 x 160) are shown. In the basis which diagonalises the Hamiltonian, the chosen 
submatrix contains elements corresponding to the dressed states closest to the cavity resonance. 



79 and 82 are the two states from the second manifold 
(two-photon excitations), closest to the resonance. At 
the positions 78 and 83 are two remaining states from the 
first manifold; at the positions 77 and 84 are two states 
from the third manifold (three-photon excitations), clos- 
est to the resonance, and so on. The states we have just 
identified suffice to indicate the dynamics of the system. 
Off-diagonal elements are coherences between the appro- 
priate dressed states. 

One striking feature can be noted immediately from 
Fig. ^. Namely, the first 'square' of elements (4 elements 
on each side) encircling the centre square (Stark states, 
2x2) remains much smaller than all of the other acces- 
sible states. In terms of dressed states, this means that 
the second manifold states remain unpopulated and the 
coherences between these and all the other states van- 
ish as well. Furthermore, as the driving increases, the 
states whose coherences with the Stark states increase 
to a noticeable size are the third manifold states (see 
Fig. | for S p = 0.09k). As the driving increases further, 
the remaining two first manifold states (their populations 
and coherences with the Stark states) also start to con- 
tribute. So, contrary to expectations, the contribution 
to the dynamics of the states beyond the effective two- 



level system do not increase according to the smallness 
of their detunings from the bare cavity frequency. The 
two least-detuned second-manifold states are in fact es- 
sentially unpopulated. 

From the above analysis of the density matrix ele- 
ments, we can deduce the solution to the 'photon statis- 
tics puzzle' of Figs. [2] and The strong photon anti- 
bunching at small £ p is the consequence of the effective 
two-level system, and has been well understood. The 
sharp rise in g^ 2 \0) with increasing E v can be attributed 
to the two-photon transitions needed to populate the 
third manifold states from the first manifold, and then 
decay back in cascade to the first manifold. These two- 
photon decays cause the sharp increase in the second or- 
der output field correlations. But, what causes the sys- 
tem to skip second manifold states? Again, the answer 
can be deduced from the density matrix. Strong coher- 
ence between the first and third manifold states, together 
with the missing population in the second manifold states 
(and vanishing coherences associated with these states) 
uncovers the effect of EIT-type quantum interference be- 
tween the dressed states. This is not surprising, since 
the quantum interference between the transitions from 
first to second manifold have been already discussed in 
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Ref. 0. 

These features are shown in Fig. ^, where the relevant 
populations and coherences are shown. Note that the 
combined population of the two inner second manifold 
states (Fig. | (b), solid line) vanishes for a large interval 
of driving, since the value of ~ 10 -16 is at the numerical 
precision value, and fluctuations are numerical, not phys- 
ical in nature. These populations become nonzero at the 
value of driving strength at which </ 2 ^(0) of Fig. || (solid 
line) peaks. These plots further justify the discussion of 
the preceeding paragraphs. 

We proceed with the development of an effective model 
with relatively few levels which nevertheless captures 
most of the important features of the dynamics. 



V. EFFECTIVE MODEL 
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In the formulation of an effective model, we rely on 
the formalism developed in Ref. |lCj. This formalism 
was very successful in the development of an effective 
two-level theory, and in explaining the fluorescence spec- 
trum. Now we extend the model and include a total of 
six dressed states in the effective model. These states are 
shown in Fig. |^. 

It should be noted, however, that the dressed states 
shown in Fig. |^ do not correspond exactly to the dressed 
states discussed in Section IV. Namely, the dressed states 



of Sec. |V| are often referred to as doubly dressed states, 
since they diagonalize the Hamiltonian with driving con- 
tributions included. In this Section, we will treat driving 
separately (see [[k| ) , and let the dressed states represent 
the eigenstates of the interaction Hamiltonian ( jla| ) alone. 
Driving can then be included through effective Rabi fre- 
quencies Vtij, coupling dressed states \d) and \ej). 

The relevant effective non-Hermitian Hamiltonian (in- 
cluding driving and damping) is thus 

= ^2P2P2 + f^3P 3 P3 + frtiPjiPjA, + ^5Pj 5 Pj5 

+ihfl i (pi -Pij 

+ik (fll 2 P2 - fi 12?4) + ifl (^13P3 - ^13^3 
+lMl 2 4 (p24 - P24) + *^25 (P25 - P25 
+iMl 3i (P34 - P34) + (P35 - P35 

-ihv p\pi - ihv 2 2p\p2 - ihT 33 p\p 3 



-ihT 44 p'l i4 p j 4 - ihT 55 p' j5 p j5 
-ihT 23 plp3 - ihT 32 p\p 2 

-ihv ib p] A p jb - ihr 54 p 



35P34 



(8) 



Operators pjk are polariton operators defined by |efc) — 
Pj k \ej)- For the states \ek) which are accessible from 
only one lower state, the notation has been abbreviated, 
so, for example, p 2 = p\ 2 . Note also that the indices 



FIG. 5: Semilogarithmic plots of the density matrix elements, 
(a) Plot of the coherences between the first manifold and third 
manifold states. Coherence is the coherence between the 
upper Stark state and lower third manifold state. Coherence 
Pig is the coherence between the lower Stark state and up- 
per third manifold state, (b) Plot of the populations of the 
dressed states. 2 nd denotes a sum of populations in the two 
relevant states of the second manifold, 1 st denotes a sum of 
populations in the two far detuned states in the first mani- 
fold and 3 rd denotes a sum of populations in the two relevant 
states of the third manifold. Parameters are as in Fig. ^, solid 
line. 



j are dummy indices, i.e. Pj^Pji — |efc)(e/|. Rabi fre- 
quencies Q12 and f2i3 have a phase term making them 
purely imaginary, while all the other Rabi frequencies are 
real. Damping terms were discussed in detail in |]Io| . It 
is now possible to formalize the distinction between the 
two dressed states bases used in this Section and Sec- 
tion IV . The basis in which the Hamiltonian (||) is writ- 
ten is defined in Ref. [flCfl and depicted in Fig. H. The 



discussion of Section IV is based upon diagonalizing the 



total Hamiltonian (including all of the 160 states used 
for the numerical simulation). 

With the effective Hamiltonian (||), we can also write 
the master equation for an effective density matrix as 



Peff = -^(fteffpeff-peffHlff^ 

+2j2SiP ef fSt, 



(9) 



1,3 



where Sk now denote the polariton collapse operators 
(Appendix |B l| ) . The effective density matrix p e f j has 
the dimension 6 x 6 - a significant reduction from 160 x 160 
used to obtain the results in Sec. IV. Equations of motion 



for the elements of p e ff are given in Appendix |B 2 . 

From the equations of motion, we can uncover terms 
which lead to the effect of quantum interference. For ex- 
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FIG. 7: Comparison of the numerical results obtained from 
the effective model including three manifolds with the results 
of numerical simulations including up to eight manifolds. De- 
tails are described in the text. The parameters are the same 
as for the solid line in Fig. ^. 



FIG. 6: Schematic depiction of the six states used to formulate 
the effective model. States |eo) and |ei) are the effective two- 
level system from Ref. juj. States |ea) and |e3) are second 
manifold states closest to the bare cavity resonance, and \e&) 
and \es) are the third manifold states closest to the resonance. 
Arrows represent effective driving. 



ample, equations for the populations of the second man- 
ifold states P22 and P33, depend on the populations ^44 
and P55 of the third manifold states and the coherences 
between these states P45 and P54; the latter with rate 
T45. At the same time, equations for P44 and P55 do not 
depend on the second manifold states, nor their mutual 
coherence. The same holds for the coherences P23 and ,045 
and their adjoints; the equation for P23 depends on P44, 
P55i P45 and P54, but not vice versa. Population of and 
coherence between the second manifold states is linked 
to the population of and coherence between the third 
manifold states. If this dependence is removed from the 
equations of motion, i.e., terms dependent on P44, ^55, 
P45 and P54 are removed from the equations for P22, P33, 
P23 and P32, cancellation of the population in the second 
manifold ceases to occur. 

The 35 equations of the effective model can, in prin- 
ciple, be solved analytically in the steady state. How- 
ever, the resulting expressions are complicated and do 
not offer significant physical insight, so we have opted to 
perform numerical solutions of the equations given in Ap- 
pendix B 2 and check for validity of the effective model. 



can be calculated as the ratio of 

(aWaa) = \w 01 \ 2 \wi 2 \ 2 P22 + k'oi| 2 |wi 3 | 2 P33 

+ [|W12| 2 |W24| 2 + |W13| 2 |W 3 4| 2 + (W* 2 -!«24 W 13W34 + C.C.)] p 4 4 
+ [|u>i 2 | 2 |w 2 5| 2 + |wi 3 | 2 |w 3 5 | 2 + (W 13 W3 5 U>12W25 + C.C.)] p 55 
+ (\w l | 2 W?12^13 P34 + C.C.) 
+ [(|W12| 2 U>241«25 + \wi3\ 2 W 3 4Wt 5 

+^12^25^13^34 + W12W24W13W35) P45 + C.C.] , (10a) 

and the square of 

(a^a) = |w i| 2 Pii + |u>i 2 | 2 p 2 2 + l w 13| 2 P33 

+ (|w 24 | 2 + |w 34 | 2 ) P44 + (\W25\ 2 + \w 3 5\ 2 ) P55 

+ [Wi2wl 3 p 2 3 + (^24^25 + w 3i w 3b) P*5 + C.CjQlOb) 

where Wij = fly/£p, and c.c. stands for complex conju- 
gate. 

The results are displayed in Fig. [?]. We have com- 
pared the numerical solutions of the effective model with 
the results of Sec. Ill and found a very good qualitative 



Having the solutions for the populations and coherences, 
the second order correlation function for zero time delay 



agreement. We do find strong bunching and threshold be- 
haviour occurring at the same values of £ p and (a^a) ss . 
However, the effective model differs from the full sim- 
ulations in the size of g^(0) at its peak, by approxi- 
mately a factor of 2. Including more states in the effec- 
tive model would yield full agreement with the numerical 
data. In particular, we found that including two states 
closest to the resonance from up to eight manifolds repro- 
duces the numerical data exactly. The reason is that the 
two-photon cascade decay responsible for the behaviour 
of g( 2 \0) can result from the decay of the higher states 
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to third manifold states first, i.e. two-photon cascade 
can, in the manner of speaking, be driven "from below" 
and "from above" . Naturally, the decay of higher lying 
states introduces more single-photon transitions as well. 
Therefore, adding one manifold at a time to an effec- 
tive model reveals that the increase in maximum value of 
g( 2 )(0) gradually diminishes with new manifolds added, 
settling at its maximum value after the inclusion of the 
eighth manifold. On the other end, the effective model 
of Fig. ^ is the smallest possible model which (at least 
qualitatively) reproduces the strong bunching effect in 
this system. 



VI. DYNAMICS OF THE FORWARD 
SCATTERING OF LIGHT 

To obtain a different and useful perspective on the 
physical processes involved in the changing nature of the 
statistics of light emitted by the coupled atom-cavity sys- 
tem, we can split the field operator into contributions 
from a coherent mean amplitude and from an incoherent 
part @, 
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Aa. 



(11a) 



where a = (a) denotes the coherent amplitude compo- 
nent of the intracavity field, while Aa denotes the inco- 
herent amplitude component (Aa) = 0, whose emergence 
is the result of scattering of the cavity field by the atom. 
Using this decomposition, one can rewrite the expression 
for (0) as 



-,(2) 



(0)-l 



(a* a* aa) ss 

(: (a*Aa + oAat) 2 ;)„ 
(M 2 + (AatAa) ss ) 2 
4|a|Re [(A a t 2 Aa) ss ] 
(|a| 2 + (AatAa) ss ) 2 



+ 



+ 



(Aat 2 Ac 



((AatAc 



(H 2 + (AatAa) ss r 



= S(Aa) + T(Aa) + V(Aa) , (lib) 

where :: denotes normal ordering, and (a^a) ss = \a\ 2 + 
(Aa^Aa) ss . The three terms in this expansion, denoted 
by S, T and V, have been identified by Carmichael 10 
for the case of a two-level atom. The decomposition (|llf ) 
shows how the behaviour of g^(0) for different values 
of driving field can be interpreted as the effect of self- 
homodyning between the coherent and incoherent com- 
ponents of the intracavity field |l{J]. From this view- 
point, it is easy to identify S(Aa) as a term describing 
the squeezing in the field quadrature in phase with the 
driving field, V(Aa) gives the variance in the incoher- 
ent component, and T(Ao) describes intensity-amplitude 
correlations in the incoherent component. Both V(Aa) 
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FIG. 8: Contributions to the correlation function g^ 2 \0) from 
the incoherent component of the intracavity field. Parameters 
are the same as in Fig. H 



and T(Aa) are determined by the correlations in the 
intensity, so the departure from the coherent value of 
the correlation function can be assigned to the effects of 
squeezing and the effects of intensity correlations. 

The contributions from the incoherent component of 
the intracavity field are shown in Fig. ^. The squeezing 
and intensity correlation parts are shown separately. The 
antibunching for weak fields comes from the squeezed 
fluctuations, which reduce the contribution from the co- 
herent scattering. As the driving increases, the squeez- 
ing decreases, but the variance in intensity fluctuations 
becomes negative, so the remaining antibunching comes 
from the sub-Poissonian intensity fluctuations in the in- 
coherent component of the field. For £ p ~ k, the 
squeezing contribution effectively vanishes, while T and 
V become positive, and antibunching disappears. Strong 
bunching clearly originates in the super-Poissonian in- 
tensity correlations, and the correlation function is dom- 
inated by V. 



VII. NONCLASSICAL BEHAVIOUR OF THE 
CORRELATION FUNCTION 

The presence of nonclassical effects in a driven atom- 
cavity system is a topic that has received much attention 
for many years. The system usually studied has been the 
canonical system of quantum optics - a single two-level 
atom confined in an optical cavity. Photon statistics in 
the bad-cavity limit was thoroughly studied by Rice and 
Carmichael [^0| , who analyzed the sub-Poissonian statis- 
tics and photon antibunching in the cavity transmission, 
for the case of weak driving. Their analysis was extended 



9 



by Carmichael et al. |EE to a system containing N two- 
level atoms. This analysis was further refined by Brccha 
et al. 0] . Clemens and Rice |^3j have extended the con- 
sideration involving a single atom to include arbitrary 
driving field strength and dephasing. In their analysis, 
Clemens and Rice pay special attention to nonclassical 
effects known as 'undershoots' and 'overshoots'. These 
are related to the violation of inequalities that hold for 
classical correlations, in particular violations that occur 
not in the value of g^(r = 0), but for certain time de- 
lays r > 0. The explanation for the undershoots has 
been given by Carmichael et al. pH in terms of quantum 
interference of probability amplitudes and collapse of the 
wavefunction. 

The requirements for the classicality of the field cor- 
relations can be derived from the Cauchy- Schwartz in- 
equality (see p3| and references therein), and expressed 
in terms of the second-order correlation function as 



£ p /k = 0.95 



£ p /k = 1.05 



|5 (2) (r)-l|<| 5 (2) (0)-l| 



(12) 



Values in excess of those allowed classically are called 
overshoots, while values below are called undershoots. 
Overshoots have been observed recently by Mielke et 
al. @. 

In this context, it is of interest to see if the over- 
shoots and/or undershoots can be found in the single- 
atom EIT-Kerr system under consideration. Photon an- 
tibunching, as an example of nonclassical photon statis- 
tics, has already been predicted H, and the effective two- 
level behaviour analyzed (|, pToj, [12|], In the present arti- 
cle, we have shown how the effects of self-homodyning of 
squeezed dipole radiation yields photon antibunching in 
the low to moderate driving limit. We have also shown 
that quantum interference between the probability ampli- 
tudes contributes to both strong antibunching and strong 
bunching, for weak and strong driving fields, respectively. 
Given this range of behaviours, we might therefore expect 
undershoots and overshoots to also occur in the single- 
atom EIT-Kerr system under suitable conditions. 

Fig. [| shows correlation functions for several values 
of driving field strength. The values of driving have 
been chosen where nonclassical behaviour is expected 
to be found. For weak driving, where the antibunching 
is strong, the delay-time dependence of the correlation 
function is well-understood in terms of the effective two- 
level system. The interesting region is for those values of 
driving for which g^ 2 '(0) increases through 1, the value 
for a coherent field. This is also the region in which the 
dynamics is well-described in terms of quantum interfer- 
ence and increased incidence of two-photon emissions. 

Not surprisingly, this is also the range of parameters 
where largest violations of the classical inequalities oc- 
cur. Since the increase in the value of g^ 2 \0) is due to 
the purely quantum effect of interference between prob- 
ability amplitudes, non-classical behaviour of the corre- 
lations can be expected. The explanation of these effects 
given by Rice and Carmichael [^0| and Carmichael et 
al. (^l), although in a different context, still holds. As 




FIG. 9: Second order correlation functions vs. time delay 
for a different values of driving field strength. Shaded areas 
denote classically allowed values, calculated from Eq. (|l2|). 



shown in Sec. VI, self-homodyning of squeezed dipole ra- 



diation with the driving field occurs in the EIT-Kerr sys- 
tem in a similar manner to that for a two-level atom. An 
alternative explanation in terms of quantum interference 
of the driving field with the atomic polarization after the 
collapse of the wave function upon a photon detection 
event offers even more insight. This is best understood 
in the context of quantum trajectory theory. 

This explanation is illustrated in Fig. [l^, where single 
trajectory realizations for the intracavity field and pho- 
ton number are plotted. We see that the collapses tend to 
occur in succession before the system returns to a (quasi-) 
steady state (which for this region of parameters happens 
after time ij 1 )- Once it returns into a steady state, a few 
Rabi cycles pass before the next set of collapses occur. 
The value of g^ (0) is determined by the ratio of the num- 
ber of jumps upwards to number of jumps downwards in 
photon number, where jump upwards suggests that the 
detection of a photon increases the probability of detect- 
ing a second photon immediately afterwards. Naturally, 
at the value of driving where g^(0) peaks (see Fig. |^), 
collapses are almost exclusively upwards, as illustrated in 
Figure O. Undershoots appear as the consequence of a 
change in sign that the amplitude undergoes at the col- 
lapse |2l| . As the system returns towards its steady state 
the polarization becomes close in magnitude and oppo- 
site in sign to the driving field, producing a near-zero 
mean intracavity field, leading to the reduced detection 
probability for a second photon. 

The overshoots can also be explained in terms of the 
collapses of the wave function. The detection of the 
first photon, emitted from the steady state situation, col- 
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small, since it is proportional to the mean intracavity 
photon number. For stronger driving fields the proba- 
bility for subsequent collapses increases, specifically due 
to the large correlations between first and third manifold 
states, as shown in Fig. ||. Therefore, the probability 
for a second photon detection after some time r < j^ 1 
increases as well, causing the correlation function over- 
shoot. In experiment, such event pairs are the source 
of delayed coincidence counts. A third collapse is also 
likely to occur before the steady state is reached. How- 
ever, the overshoot disappears (or significantly decreases) 
for all time delays r after the first peak. This is expected, 
since the exact form of the wave function after the second 
collapse depends on the delay time between the second 
and third photon. Averaging over all possible realiza- 
tions washes out the nonclassical effects due to different 
possible evolutions following the second collapse. 

A stronger driving field causes more subsequent col- 
lapses to occur, and nonclassical correlations are com- 
pletely washed out. We find that overshoots and under- 
shoots vanish at driving strength £ p 1.14k. 



VIII. CONCLUSION 



FIG. 10: Intracavity photon number and field amplitude in a 
typical realization of a single quantum trajectory, for £ p /k = 
l. 




time 



We have presented an analysis of the properties of the 
photon statistics of light emitted by a single atom intra- 
cavity EIT-Kerr system. It was found that the statistics 
change abruptly as the driving field strength increases. 
Specifically, strong photon antibunching, dominant in the 
weak to moderate driving regimes, is replaced by a strong 
photon bunching in the output field for the stronger driv- 
ing. We have identified the effect of quantum interference 
between the dressed states to be responsible for this sud- 
den change, and presented an effective model explaining 
qualitative features of this behaviour. 

Furthermore, we have analyzed contributions from the 
incoherent scattering to the system dynamics and found 
that the strong photon antibunching can be explained 
in terms of the self-homodyning of the incoherent intra- 
cavity component with the coherent component for weak 
driving, and in terms of reduced intensity fluctuations 
for moderate driving strengths. Strong bunching is the 
signature of super-Poissonian intensity fluctuations. 

Finally, nonclassical behaviour of the correlation func- 
tion was found, and the effects of undershoots and over- 
shoots analyzed. 



FIG. II: Intracavity photon number in a typical realization 
of a single quantum trajectory, for £ p /k = 1.75. 
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APPENDIX A: STARK SPLITTING OF THE 
DRESSED STATES 



can write field and atomic collapse operators in terms of 
the bare states as 



In this Appendix, we briefly review the effect of dy- 
namic Stark splitting of dressed states, fully elaborated 



upon in Ref. 10 . 

Dressed states of the Hamiltonian (|la|) contain, among 
others, two states on exact cavity resonance, separated by 
energy hiu cav . These are the ground state and a state be- 
longing to the first excitation manifold. In the schematic 
depiction of Fig. |], these are denoted by |eo) and |ei). 
Driving by the external field £ v enters through the ef- 
fective Rabi frequency Qqi — £ p / \A + g\l^ 2 c - Resonant 
coupling of the two levels causes splitting of the two states 
into a symmetric and antisymmetric linear combination 
of the two 



|V±) = (|eo>±|ei))/V2, 
with corresponding energies 



e± = ±Vng x - (r /2)' 



(Al) 



(A2) 



where T = k/ (l + gf/ft^j is the decay rate of the ex- 
cited state |ei). The splitting of energy levels occurs at 
the threshold values of driving field 



k/2 



VT+gJM 



(A3) 



The two states \ip±) are the two states closest to the 
cavity resonance for the range of driving strengths con- 
sidered in this article. 



APPENDIX B: QUANTUM JUMPS AND 
DENSITY MATRIX ELEMENTS 

In this Appendix, we derive jump terms for the effec- 
tive master equation (^) , and the equations of motion for 
density matrix elements. 



Cl2 

C32 
034 



|0, 1)(1, 1| + |0, 3)(1, 3J -r- jO, 4)(1, 4| 
W2(|l, 1>(2, 1| + |1, 2) (2, 2| + |1, 3)(2, 3|) 
+V3|2, 1><3, 1] , (B2a) 
|1, 1)(1, 2| + |2, 1)(2, 2|, (B2b) 
|1, 3)(1, 2| + |2, 3)(2, 2|, (B2c) 
|0, 3)(0,4| + |1, 3)(1, 4|. (B2d) 



It can be deduced from the bare states that 
the collapse operators can be expressed as lin- 
ear combinations of the following polariton op- 
erators: a cx {pi, p 2 , P3, P24, P3A, P25, P35}; 
CT12, cr 34 cx {p 2 , P3, P24, P34, P25, P35} and a 32 oc 

[p 2 P2, P\P31 PiPi, PIP5, P 2 P3, p\p2, Pj4Pj5, p)bPj±)- 

Note that the polariton jump terms arising from the 
atomic jumps associated with the operator (732 (i.e. 
proportional to 72) couple dressed states within the 
same manifold. This is possible since the associated 
atomic transition is not coupled to the cavity mode, but 
to the classical field il c , and therefore jumps occurring 
in this atomic decay channel do not switch between the 
adjacent manifolds. 



Using the correspondences above, one can write all of 
the jump terms in the polariton basis, thus making the 
transition Sp — ► iS e //p e //. This transformation can be 
viewed as a rotation of a truncated basis of a superoper- 
ator. The coefficients Tijki of the jump terms are given 
in Table | 



1. Jump Terms 



We follow the notation of Carmichael |T^] and rewrite 
the effective master equation (^) in a Lindblad form as 
p e ff = C e ffp e ff, where the Liouvillian superoperator 
can be divided into a part describing the free evolution 
between the jumps, C e ff — Seff, and a part describing 
the jumps, S e ff. Jump term S e ffp e ff arises from the 
equivalent jump term in the full master equation, 



Sp = 2k apa 1 ^ + 271 (T12PCT21 

+272 CT32PC23 + 273 034/0043 



(Bl) 



Using the methods of Ref. |l(j] , and consistent with the 
truncation of the dressed states space (as in Fig. ||), we 



2. Equations of Motion 



The equations of motion for the density matrix ele- 
ments in the basis schematically shown in Fig. || can be 
found from the master equation (||). There is a total of 
35 equations, since population conservation can be used 
to eliminate one of them. The equations for populations 
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P03 



P04 



P05 = 



Pl2 



POO 

Pn 



P22 



P33 



Pa 



P55 



2r pu + Qo (poi + P10) , 

— 2TqPu + Tii22P22 + Tn33P33 
-Qq (poi + pio) 

+f2l 2 pi2 + ^12/921 + ^13Pl3 + ^*3P31 , 

— 2 (T22 + T2222) P22 

-r 2 3 (P32 + P23) + 2r 45 (p 45 + p 54 ) 

— Ol2Pl2 — ^*2P21 + T2233P33 + I^2244P44 + 
+ ^24 (P24 + P42) + ^25 (P25 + P52) , 

— 2 (T33 + T3333) P33 

-r 2 3 (^32 + P23) + 2r 45 (p 45 + p 54 ) 

— ^13Pl3 — ^13P31 + T3322P22 + ^344/544 + 
+ ^34 (P34 + P43) + ^35 (P35 + Pbz) , 

— 2 (T44 + T4444) P44 + T4455P55 — 2F45 (^45 

-fi 2 4 (P24 + P42) - ^34 (P34 + P43) ) 

— 2 (r 55 + r 5 55 5 ) p 55 + r 55 4 4 p44 — 2r 45 (p 45 

-^25 (P25 + P52) - ^35 (P35 + P53) ■ 



(B3a) 



(B3b) 



Pl3 



r2255P55 P14 

(B3c) 

Pl5 

r3355P55 p23 

(B3d) 

+ P54) 

(B3e) 

+ P54) p24 

(B3f) 

P25 



P34 



The equations for coherences are 



— (^33 + r 33 33 — 163) po3 — T 2 3P02 
+ro314Pl4 + ro3i5pi5 

+f2 pi3 - ^13P01 + ^34P04 + ^35P05 (B4c) 

— (T44 + T4444 + Z£ 4 ) p 4 — T45P05 

~^24P02 ~ ^34P03 + ^0Pl4 (B4d) 

— (T55 + r 5 555 + ie 5 ) po5 - r 45 p 4 

~^25P02 ~ ^35P03 + ^0Pl2 (B4e) 

— (To + T22 + T2222 — fa) Pl2 — T23P13 
+ri224P24 + Ti225p25 + Ti234p34 + Ti235P35 
+ ^12 (P22 - Pll) 

-f2 P02 + ^i 3 P32 + ^24Pl4 + ^25Pl5 (B4f) 

— (To + T33 + r 3 333 - ie 3 ) P13 - T 2 3Pi2 
+ri324P24 + Ti325P25 + ri334P34 + 1^1335/535 
+ ^13 (P33 - Pll) 

-f\)P03 + ^l 2 P23 + ^34Pl4 + ^35Pl5 (B4g) 

— + T 4 4 + T4444 — 164) pu — r 45 pi 5 — r2oP04 
+0 2 4Pl2 + ^34Pl3 + ^i 2 P24 + ^ 3 P34 (B4h) 

— (r^O + ^55 + T5555 — ie$) P15 — T45pi4 — £1()P05 
+f2 2 5Pl2 + ^35Pl3 + ^l 2 P25 + ^ 3 P35 (B4i) 

— [T22 + T2222 + T33 + 1^3333 + T2233 — i fa + £3)] P23 

— 1^23 (P22 + P33) + T2344P44 + T2355P55 + T2332P32 
+r2345P45 + T2354P54 + f^ 2 Pl3 + ^*3P21 

+^34P24 + ^35P25 + ^24P43 + ^25P53 (B4j) 

— [r 2 2 + r 2 222 + r 44 + r4 444 + r 2 244 — i fa + £4)] P24 

— T23P34 — T45P25 + T2435P35 

-^24 (P44 - P22) - ^12Pl4 - ^34P23 + ^25P54 (B4k) 

— P22 + T 2 222 + r 55 + r 55 55 + r 22 55 - i fa + £5)] P25 

— T23P35 — T45P24 + T2534P34 

-^25 (P55 - P22) - ^12Pl5 - ^35P23 + ^24P45 (B41) 

; (£3 + £4)] P34 



r 



3344 



P35 



P45 = 



— T33 + T3333 + T44 + T4444 

— T23P24 — T45P35 + T3425P25 

-^34 (P44 - P33) - ^13Pl4 + ^35P54 ~ ^24P32 (B4m) 

— ^33 + T3333 + T55 + 1^5555 + 1^3355 — i fa + £5)] P35 

— T23P25 — T45P34 + T3524P24 

-^35 (P55 ~ P33) ~ ^13Pl5 + ^34P45 ~ ^25P32 (B4n) 

— T44 + r 444 4 + r 55 + r 5 555 + r 4 4 55 — ifa + £ 5 )] p4 5 

~^45 (P44 + P55) + T4554P55 

~^24P25 ~ ^34P35 ~ ^25P42 ~ ^35P43 ■ (B4o) 



Damping coefhcients To and have been evaluated in 
general in Ref. (To). Coefficients Tijki are related to the 
jump operators and are given in Table ffl. The source of 
x _ p „ 1 op „ 1 or „ the different terms in these equations should be apparent 

POI — — -L 0P01 + *L 02P12 + ^ 03P13 . . ,. .. . -rC , 1 . 

^ . in ft fB4 ) horn the earlier discussion, h or example, damping emerg- 

+ (P11 Poo J + 12P02 + 13P03 ( a J ing from jump terms is denoted as Tij-fe; , while damping 
P02 = — (r22 + T 2 222 — fa) P02 — T23P03 coming from the effective Hamiltonian (||) is denoted as 

+r 214Pl4 + r 02 15Pl5 To and Tij. 

+^oPi2 — ^*2Poi + ^24Po4 + ^25Po5 (B4b) It is straightforward to prove that the following equal- 
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ities hold: 

2(r 2 2+r 2 222) = rn22+r 33 22, (B5a) 

2 (r33 + ^3333) = r 1133 + r 2 23 3 , (B5b) 

2 (F44 + F4444) = T 2 244 + r 33 44 + ^544 , (B5c) 

2 (r55 + r 55 5 5 ) = r 2 255 + r 3355 + r44 55 , (B5d) 



preserving the population conservation condition 
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2332 


72 ^ 2) */3« 2 V 2 2) *ft 2) 


2534 


72 Mo Po Po lit, 

I f Z 1 Z 1 O ' O 


2435 


T 2ft 2) */f ) *pf ) ft 3) 

' ' Z ' o 1 ' z 


3425 


72/4 2) */3f ) *4 2 V 2 3) 


3524 


72/4 2) *pf , *pf ) /4 3) 


2255 


n\w 2& \ 2 + 7 i|ft 3) | 3 |a 2 2) | 2 +73|M 2 2) | 2 |4 3) | 2 


3355 


^K5| 2 +7l|M3 3) | 2 |4 2) | 2 +73|ft 2) | 2 |^ 3) | 2 


2355 


* (2)*l (3)i2 (2) . (2)*i (3)i 2 (2) 
+ 7l«2 ft a, + 73ft ft 


1235 


KW13W25 + 71Q, p., a., /i., + 73ft ft za, v-i 


1234 


™ «,• +^i 1 >/9«> «»>+^ l >«PVWV« 

+ 7l a i P3 Co ft + 73ft ft ^9 

' ± ' o z ' z ' " r ' z z 


1325 


KW12W35 + ■yia^* pi 2 ^ + 73ft 1) *M3 2>J/ 3 3) * iy 2 2) 

•''J ' i 'Z O'O ''-''O 'O O Z 


2344 


* , (2)* (3) (2) (3)* . (2). (2), (3), 2 
KIII24H34 + 7l Q 2 ft Q 3 ft + 73ft ft ^2 

ii^ ' Z ' Z O ' Z ' ^ / z ' O I Z 1 


2345 


* , (2)* (3) (2) (3). . (2)* (2) (3) (3)* 
KW 2 iW 35 + 7ia 2 M2 Q 1 ft + 73ft ft v 1 v \ 

-* 1 -*- ^ »Z O'O ; v ^z 'O Z O 


2354 


* , (2)* (3)* (2) (3) . (2)* (2) (3)* (3) 

KW25W34 + 7ia 2 u 2 a, /ii + 73ft ft VA ^3 

' Z ' Z O'O f v /z 'O Z O 


4554 


(3)*a(3).a(3) (3) 
72Mo P9 P3 ft 

1 ' Z ' Z 1 O ' 


2222 


( 7 2/2)| M2 2, | 2 |/3 2 2, | 2 


3333 


(72/2)!Mf| 2 |4 2, | 2 


4444 


( 7 2/2)| M2 3, | 2 |/3f , | 2 


5555 


(72/2)! M < 3) | 2 |pf | 2 


2233 


( 72 /2)| M2 2, | 2 |/3f | 2 


2244 


.| W 24| 2 +7 1 |4 2, | 2 |ft 3, | 2 +73|ft 2, | 2 |^ , | 2 


2255 


K |™ 2 s| 2 + 7 i|a 2 2, | 2 |Mr , | 2 +73|M 2 2, | 2 |^ 3, | 2 


3344 


.K4| 2 +7i|4 2, | 2 |ft 3, | 2 +73|ft 2, | 2 |4 3, | 2 


3355 


K |™3 5 | 2 +7i|4 2, | 2 |ft 3, | 2 +73|ft 2, | 2 kf , | 2 


4455 


( 72 /2)| M2 3) | 2 |/3< 3) | 2 



TABLE I: Expressions for the 'jump' coefficients. Here, Wij = £lij/£ p , and the parameters a, [3, /i and v are given in lid 



